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J2lc Abstract 

[ The ground state properties of interacting Bose gases in external 

^ I potentials, as considered in recent experiments, are usually described 

M ■ by means of the Gross-Pitaevskii energy functional. We present here 

^ , the first proof of the asymptotic exactness of this approximation for 

I the ground state energy and particle density of a dilute Bose gas with 

^ • a positive interaction. 

1 Introduction 



Recent experimental breakthroughs in the treatment of dilute Bose gases 
have renewed interest in formulas for the ground state and its energy derived 
many decades ago. One of these is the Gross-Pitaevskii (GP) formula for the 
energy in a trap ||l]]-0, such as is used in the actual experiments. We refer to 
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1^ for an up to date review of this approximation and its applications. One 
of the inputs needed for its justification is the ground state energy per unit 
volume of a dilute, thermodynamically infinite, homogeneous gas. This latter 
quantity has been known for many years, but it was only very recently that it 
was derived rigorously P for suitable interparticle potentials. Consequently, 
it is appropriate now to use this new result to go one step further and derive 
the GP formula rigorously. 

The starting point for our investigation is the Hamiltonian for identical 
bosons 

N 

H(N) = ^ (_v2 + V{x,)) +Y,v{\x, - (1.1) 

i=l i<i 

acting on totally symmetric, square integrable wave functions of (xi, . . . , xat) 
with Xi G M^. Units have here been chosen so that fi = 2m = 1, where m is 
the mass. We consider external potentials V that are measurable and locally 
bounded and tend to infinity for — > oo in the sense that mi\x\>R V{x) oo 
for R oo. The potential is then bounded below and for convenience we 
assume that its minimum value is zero. The ground state of — + V{x) pro- 
vides a natural energy unit, hw, and the corresponding length unit, ^Jhjmuj, 
describes the extension of the potential. We shall measure all energies and 
lengths in these units. In the available experiments V is typically ~ |a;p and 
\Jh/mio of the order 10~^ m. 

The particle interaction v is assumed to be positive, spherically symmetric 
and decay faster than at infinity. In particular, the scattering length, 
denoted by a, should be finite. We recall that the (two-body) scattering 
length is defined by means of the solution u{r) of the zero energy scattering 
equation 

-u"{r) + \v{r)u{r) = (1.2) 

with m(0) = 0; by definition, a = limr^oo(?' — u{r) / u' {r)) . Let fi(r) be 
a fixed potential with scattering length ai. Then v{r) = {ai/a)'^vi{air /a) 
has scattering length a. We regard in the following vi as fixed, but vary a 
(in fact, we shall take a = ai/N). The ground state energy E'^^ of ( |1.1D 
depends on the potentials V and w, besides N , but with V fixed and v{r) = 
{ai/a)'^vi{air/a), the notation E^^{N,a) is justified. The corresponding 
eigenfunction will be denoted ^I/q^'*. It is unique up to a phase that can be 
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chosen such that the wave function is strictly positive where the interaction 
is finite 0. The particle density is defined by 

pTai^)=N [ \¥^\x,X2,...,XN)\'d'x2...d'xN. (1-3) 



The Gross-Pitaevskii (GP) energy functional is defined as 

£GP[$] = /■ (|V$(x)|2 + V^(x)|$(x)|2 + 47ra|$(x)|^)c/3x (1.4) 

where $ is a function on M^. For a given the corresponding GP energy, 
denoted E^^{N, a), is defined as the infimum of £^[<l>] under the normalization 
condition 

[ \<^{x)\'^d^x = N. (1.5) 

It has the simple scaling property 

E^^{N, a) = NE^^{1, Na). (1.6) 



What ( |1.6| ) shows is that the GP functional ( |1.4| ) together with the normal- 
ization condition (|1.5|) has one characteristic parameter, namely Na. (Recall 
that lengths are measured in the unit ^^h/muj associated with \^ so a is 
dimensionless.) Thus, if we want to investigate the non-trivial aspects of GP 
theory we have to consider a limit in which N oo with Na fixed. This 
explains the seemingly peculiar limit in Theorems 1.1 and 1.2. As Na — »• oo 
the GP energy functional simplifies, since the gradient term becomes small 
compared to the other terms, and the so called "Thomas-Fermi limit" de- 
scribed in Theorem p.2| results. In some typical experiments a is about 10"'^, 
while varies from 10^ to 10^ Thus Oi in Theorems 1.1 and 1.2 varies from 
1 to about 10^. 

In the next section it will be shown that the infimum of the energy func- 
tional (|1.4|), under the subsidiary condition (|1.5|), is obtained for a unique, 
strictly positive function, denoted The GP density is given by 
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It satisfies 

pZ(^) = NpfU^). (1.8) 

The main result of tliis paper concerns tlie befiavior of the quantum me- 
chanical ground state energy E^^{N,a) when N is large, but a is small, so 
that Na is 0(1). It is important to note that although the density tends to 
infinity for ^ oo (by Eq. (|1.8|)) we are still concerned with dilute systems 



in the sense that a p ^ 1, where 

P=^l pTMd'x (1.9) 
is the mean GP density. (Note the exponent 2 in ( |1.9| ).) In fact, since 

The precise statement of the limit theorem for the energy is as follows. 

Theorem 1.1 (The GP energy is the dilute limit of the QM energy). 

For every fixed ai 

...n ^ ^-(MO (1,10) 

N^oo iV 

and the convergence is uniform on bounded intervals of ai . 

While we do not prove anything about Bose- Einstein condensation, which 
necessarily involves the full one-body density matrix p^^\x,x'), we can make 
an assertion about the diagonal part of the density matrix, p*^^(x) = 
p^^^x, x): 

Theorem 1.2 (The GP density is the dilute limit of the QM density). 

For every fixed ai 

in the sense of weak convergence in L^. 



For the proof of Theorem |TTT| we establish upper and lower bounds on 
E^^{N, a) in terms of E'^^ {N, a) with controlled errors. Theorem |1.2| follows 



from Theorem |1.1| by variation of the external potential. The upper bound is 
obtained in Section 3 by a variational calculation which generalizes the upper 
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bound of Dyson for a homogeneous gas of hard spheres. We also derive 
an upper bound on the chemical potential, i.e., the energy increase when one 
particle is added to the system. This upper bound is used in the proof of the 
lower bound of the energy in Section 4. The main ingredient for the lower 
bound, however, is the bound for the homogeneous case established in p. 
In addition, some basic properties of the minimizer of the GP functional are 
used in the proof and we consider them next. 

2 The Gross-Pitaevskii Energy Functional 

The GP functional is defined by (|1.4|) for $ G I' with 

D = {$ : V$ G L2(M^), V"|$|2 G L\W^), $ G L^(M^) H ^^(M^)}, (2.1) 

where / G /.^(R"") means J^„ < oo. The corresponding GP energy 

is given by 



E^P(Ar, a) = inf{£^P[<l>] : $ G Pjv} (2.2) 



with 



VN = Vn{^: J\^{x)\'^d^x = N}. (2.3) 



The basic facts about the GP functional are summarized in the following 
theorem. 

Theorem 2.1 (Existence and properties of a minimizer). The infi- 
mum in (|2.2| ) is a minimum, i.e., there is a G "Pat such that 

E'^^{N,a) = S^^[^'^^]. This is unique up to a phase factor, which can 
be chosen so that ^'^^ is strictly positive. is at least once continuously 
differentiable, and if V is C°° then also $*^p is C°° . The energy i?^P(A^, a) 



is continuously differentiable in a and hence (by Eq. ( \l.di) ) also in N. The 
minimizer solves the Gross-Pitaevskii equation 

-V^$ + 1/$ + 87ra|$p$ = /i$ (2.4) 

(in the sense of distributions) with 

/i = dE^^{N, a)/dN = E^^{N, a)/N + A-nap. (2.5) 

Here p is the mean density ( |1.9|) . 
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The GP energy functional is mathematically quite similar to the energy 
functional of Thomas-Fermi-von Weizsacker theory and Theorem |2.1| can be 
proved by the methods of Sect. VII in p. For completeness, the proof is 
given in Appendix A. With additional properties of V one can draw further 
conclusions about 

Proposition 2.1 (Symmetry and monotonicity). // V is spherically 
symmetric and monotone increasing, then ^"^^ is spherically symmetric and 
monotone decreasing. 

Proof. Let $* be the symmetric-decreasing rearrangement of (see 
Then£GP[$*] <^GP[$GP]_ □ 



Proposition 2.2 (Log concavity). // V is convex, then is log con- 
cave, i.e. $*^P(x)^<l>'^^(y)(^"'^^ < (Ax + (1 - A)y), for all x,y G M^, 
AG (0,1). 

Proof. Using the Trotter product formula it suffices to show that the solutions 
u{t, x) of the equations 

— V M = 0, — + Vu = 0, — + Snau^ = uu 

ot at ot 

are log concave, if m(0,x) is log concave. The first follows from the fact, 
that the convolution of two log concave functions is log concave, the second 



follows easily from convexity of V^, and the third is shown in □ 



The GP theory has a well defined limit if Na oo. It is sometimes 
referred to as the "Thomas-Fermi limit" of GP theory because the gradient 
term vanishes in this limit. For simplicity we restrict ourselves to homoge- 
neous external potentials V, i.e., 

V{Xx) = X'V{x) (2.6) 

for some s > 0. 

Theorem 2.2 (Large Na limit). Let V be homogeneous of order s and let 
T he the functional 

T\p\ = I [V{x)p{x) + ATiap{xf) d^x (2.7) 
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with p{x) > 0, X & R^. Let F{N,a) be the infimum of T under the 
condition J p = N. By scaling, F{N,a) = NF(l,Na) and F{l,Na) = 
(7Va)*/(*+^)F(l, 1). In the limit Na ^ oo we have 



(2.8) 



The minimizing density of T under the condition J p — 1 and with a — 1 is 
given by 



pUx)^i8nr'[p-Vix)]_ 



(2.9) 



with p = F(l,l) + 47r J (pi^i)^, and = t for t > and otherwise. 
Moreover, 



Na—^oo 



(2.10) 



strongly in L? 

Proof. Since S^^[^] > J^[p] it is clear that E^^{l,Na) > F{l,Na). For 
the converse we write p in the form p{x) — (A'"a)~^/^*+^^p((A^a)~^/(*+^^x) 
and obtain 



+ Vp + Attp^ ) d-^x 



JP-[p] = (A^a)^/(^+^) J {Vp + ATrp')d^x. 
In particular, F{l,Na) — (A^a)*/(^+^)F(l, 1), and with p — we obtain 



< oo.) In the 



(Regularizing V, if necessary, we may assume that J Vypf 

limit Na — > oo the gradient term vanishes, and thus the limit of the energies 
is proved. Now 



lim , X // o^ 

iP 1 

I - ir i_,Na\ 

hm 



{Na} 



a + 2 
s+3 



F(l,l). 
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Since 1) is the minimum of T / {Nay^^^^'^\ the second term vanishes 
for Na — > cxD, and it follows that pf^^ is a minimizing sequence for 
J {y p + 4:71 p^) . Since both terms in the functional are nonnegative, they 
must converge individually, in particular HpfAralh converges to ||pf ilh- On 
the other hand p^^^^ converges weakly to by uniqueness of the minimizer. 
Together with the convergence of the norms this implies strong convergence. 
The solution of the variational equation for pf ^ is simply pf i = (87r)~^[p— 

V]+ with p given by p = 1) + 47r / (pf,i)^ □ 



Lemma 2.1 (Virial theorem). When V is homogeneous of order s, as in 



( \2. dj ), the minimizer of the GP functional satisfies 

^ j |V$^P(x)prf3^-| j $^P(x)V(x)rf3x + 47ra j ^^^(x)^^^^ = 0. 

(2.11) 

Proof. Define $fc by 

Because is the minimizer of it must be true that 

d 



= 0. 

fc=i 



dk 

This leads to the virial theorem ( p.ll|) . □ 



In the proof of the lower bound we shall also consider the GP energy 
functional in a finite box. For i? > we denote by A/j a cube centered at the 
origin, with side length 2R. The energy functional E'g^ in the box is simply 
(|1.4| ) with the integration reduced to A/j, the corresponding minimizer, de- 
noted by satisfies Neumann conditions at the boundary of A/;, and is 
strictly positive. Analogously to (|1.7|) , (|1.9|) and (|2.5| ) we define p^^, p/j and 



Pr. The corresponding energy will be denoted by -E^^(A^, a). Then we have 
the following 

Lemma 2.2 (GP energy in a box). 

hm E^P(iV,a) = E^P(Ar,a) (2.12) 
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Proof. Using N^^'^^'^^ xr/ W^'^^ Xr\\2 as a test function for , where xr 
denotes the characteristic function of A/j, we immediately get 

hm Ef {N, a) < E^^ {N, a). (2.13) 

R^oo 

Let Qr be a C°° function on M^, with Qu = outside Ar, and 6^ = 1 inside 
A/j-i. We use N^^'^^fOR/W^fQRh as a test function for . Since VGr 
is bounded and 

hm / \^7\' = 

R^oo J 

Ar\Ar-i 



(because V tends to infinity and ] is bounded by (|2.13| )), we have 



E^^{N,a) < hminf £;^P(Ar,a). 

R^oo 



This completes the proof of (|2.12|) . □ 



3 Upper Bounds 

3.1 Upper bound for the QM energy 

It will now be shown that for all and small values of ap^^^ (with p = 
^ p'^^{xfd^x/N, cf. (|r|)) 

EQ^(Ar, a) < E«P(Ar, a)(l + (ap^)) . (3.1) 

This upper bound, which holds for all positive, spherically symmetric v with 
finite scattering length, is derived by means of the variational principle. We 
generalize a method of Dyson [0, who proved an upper bound for the homo- 
geneous Bose gas with hard-sphere interaction. Consider as a trial function 

^ = F{x,,...,xn)G{x,,...,xn) (3.2) 

with 

F{xi,...,xn) = llf^^Fi{xi,...,Xi), G{xi,...,xn) = U^i9ixi) (3.3) 
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where 

Fi{xi, . . . ,Xi) = f{ti), ti=mm{\xi-Xj\J = l,...,i-l), (3.4) 
with a function / satisfying 

< / < 1, /' > 0, (3.5) 

and 

^7(a;) = $^P(x)/||$^P|U. (3.6) 

The function / will be specified later. This trial function is not symmetric 
in the particle coordinates, but the expectation value (\l/|if^^)\l/)/(\E'|\E') is 
still an upper bound to the bosonic ground state energy because the Hamil- 
tonian is symmetric and its ground state wave function is positive. Hence 
the bosonic ground state energy is equal to the absolute ground state energy 

(i, 0). 

The physical meaning of the trial function can be understood as follows: 
The G part describes independent particles, each with the GP wave function. 
The F part means that the particles are inserted into the system one at 
a time, taking into account the particles previously inserted, but without 
adjusting their wave function (cf. 0]). Although a wave function of this form 
cannot describe all correlations present in the true ground state, it captures 
the leading term in the energy for dilute systems. 

For the computation of the kinetic energy we use 

[ W^^= [ {GVlG)F^- [ G2|VfcF|2, (3.7) 

where Vfc denotes the gradient with respect to x^. A: = 1, . . . , A^. We write 

{1 for i = k 
-1 for U = \xi - Xk\ ■ (3.8) 
otherwise 

Let rii be the unit vector in the direction of Xi — Xj(^i), when Xj(^i) is the nearest 
to Xi of the points (xi, . . . ,a;j_i). (Note that really depends on all the 
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points xi, . . . ,Xi and not just on the index i. Except for a set of zero measure, 
is unique.) Then 

GVkF = J2 '^F-h.knJ'iti), (3.9) 

i 

and after summation over k 



<2|*|^^i^-V'(i.)^ + 2|*p hkejk\F-'F-'f{U)f'{tj). 

i k<i<j 

The expectation value can thus be bounded as follows: 
^"=Z^ rKT/i2 



(3.10) 







-'F-' fit,) fit,) 




2 



+ 2E n,^i (3.11) 



" / l*P (-S(xi)-V?9(xi) + V{xi)) 



i=l 



For i < p, let F^^j be the value that Fp would take if the point Xi were omitted, 
i.e., 

Fp,i = f{\xp-Xk(p)\), (3.12) 

where Xk(^p) is the nearest to Xp of the points [xi, . . . ,Xi_i,Xi^i, . . . ,Xp_i). 
The reason for introducing these functions is that one wants to decouple the 
integration over Xi from the integrations over the other variables. (Note that 
Fp^i is independent of Xi.) Analogously, one defines Fp^ij by omitting Xi and 
Xj. This decouples simultaneously Xi and xj from the other variables. The 
functions Fi occur both in the numerator and the denominator so one needs 
estimates from below and above. Since 

Fp = mm{Fp^ij, f{\xp - Xj\), f(\xp - Xi\)}, (3.13) 

one has (recall that / < 1) 

Fl,jf{\xp - x,\ff{\xp - x,\f < Fl < (3.14) 
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Hence, with j < i, 

• • • ■■■Fjj < -Ff+ij . . . . . . -F^^ij (3.15) 

and 

Fj ■ ■ ■ Ffj > Fj^^j . . . F^_-^ jFI^^ ^j . . . F^ -j (3.16) 

N \ / ^ 



We now consider the first two terms in (|3.11|) . In the numerator of the first 
term for each fixed i we use the estimate 



i-l 



and in the second term we use Fi < f{\xi — Xj\). For fixed i and j one 



ehminates and Xj from the rest of the integrand by using ( 3.15 ) and 



Fj < 1 in the numerator and ( |3.16|) in the denominator to do the Xi and 



Xj integrations. With the transformation rj = Xi — Xj, X = i^i + ^i)/2 one 
gets 

(2/'(|xi - Xj\f + v{\xi - Xj\)f{\xi - Xj\f) g{xifg{xjfd^Xid^Xj = 
{2f'i\r^\f + v{\v\)fi\v\f) g{x + Wgix - \rifdSd\. 

(3.18) 

By the Cauchy-Schwarz inequality 

j 9{X + Wg{x - \vfd\ < j gix^d'x, (3.19) 

and one obtains 

{2f'{\xi - Xj\y + v{\xi - Xj\)f{\xi - Xj\f) g{xifg{xjfd^Xid^Xj 

(3.20) 

g{xtd\ / {2f\\r^\f + v{\^\)f{\^\f)d'^^2 gix^d'xJ. 



< 
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In the denominator one gets, using that < (7 < 1, 



/ (1- E (1-/(1^ 



- f{\xp - XilY) I g{xifSxi 
P=i,P^i J (3.21) 

> / gix^d'x -N [{1- f{\xp - x,\f) = [ gixfd'x - NI. 



The same factor comes from the x^-integration, the remaining factors are 
identical in numerator and denominator, and so finally the first and second 
term are bounded by 

j^^ ' {j g{xy<Px-Niy - {j g{xy<Px~Niy ^' ^ 

A similar argument is now applied to the third term of ( |3.11|) . Note that 
the contributions from k = i and k < i are the same. Therefore 



/ \e^k\f{U)f\U)\e^u\f{t,)f\tj)g{x,fg{x,f(fx4^x, < 
k=i ^ 

i — l « 

2^ / f{\xi - Xk\)f{\xi - Xk\)f{\xj - Xk\)f'{\xj - Xk\)g{xifg{xjfd^Xid^Xj. 

k=i 

(3.23) 

With g <1 one gets 

2E / fi\x^\)f\\x.\)f{\x,\)f{\x,\)d'xd' 



Xj 



k=l " 

= 2(z-l)( / f{\x\)f'i\x\)d'x] ^2{i-l)K 



2 

3^ \ — n(A 1 \ 1^2 



(3.24) 



The summation over i and j gives 



Y.Y.^^~1) = \n{N-1){N-2). (3.25) 



6 

j=i i=i 



In the denominator one gets the same factors as above, and so the third term 
is bounded by 

2 3 

3^ g{xfd^x- Nlf ^^'^^^ 
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Next consider the last term of Define e by 

{-g{x)V^g{x) + V{x)g{xf) £x = e I g{xfd^x. (3.27) 



After eliminating Xi from the integrands in the numerator and the denomi- 
nator and using Fi < 1 one sees that the term is bounded above by 

ejg{xfd^x 

jg{xyd^x-Nr ^^-^^^ 

Putting all terms together we obtain as an upper bound for the ground 
state energy 

efg(xfd^x [g(x)^d^xJ 2,. 



N - J g{xYd^x - NI (/ g{xyd^x - Niy 3 (/ g{xfd^x - Nlf 

(3.29) 

with /, J, K and e defined by [W^ . i WM . O) and (IOTI) . It remains 
to choose /. We take for h > a (we shall soon fix h) 

j^^^^[{l + e)u{r)/r for r<b ^^^^^ 
I 1 for r > b 

where u{r) is the solution of the scattering equation 

-u'{r) + ^v{r)u{r) = with u{0) = 0, lim u'{r) = 1 (3.31) 

and e is determined by requiring / to be continuous. Convexity of u gives 

1 for r < a , , , I for r < a 

r > uir) >\ - , 1 > u'ir) >\ " . (3.32) 

\r — a for r > a II — ^ for r>a 

These estimates imply 

J<{l + ef4:7ia (3.33) 
/ < 47r ("y + ab{b - a)\ (3.34) 



< 47r(l + e)a (^6 - (3.35) 

< e < (3.36) 

b — a 
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For ( |3.33| ) we used partial integration. By definition ( |1.9| ) 



2 ' 



(3.37) 



and we choose b such that 



\r II oo 



= C. 



(3.38) 



With this choice the factor in the denominators in (|3.29|) is bounded by 



NI> g^l-T 



Note that c < 1, and a < b holds provided 



a 
63 



Air 



3 II ^GP|| ^ 1 

oo ^ 



-a \\p 



(3.39) 



(3.40) 



Collecting the estimates (|3.33| )- (|3.36|) , we finally obtain 
Theorem 3.1 (Upper bound for the QM energy). 



1 + 2 /aN _ 2 /aN2 ^ r a\3 



(3.41) 



with b and c defined by ( p. 381 ). 



Remark 3.1 (Negative potentials with hard core). ( p.l|) can be ex- 
tended to include partially negative potentials of the form 



oo for < r < d 

v{r) = <^ —\w{r)\ ioT d < r < Ro 
for r > Rq, 



(3.42) 



as long as f'{rY + |t'(r)/(r)^ > for all r. With / from ( |3.30D , this is the 
case for sufficient shallow potentials. The potential energy is then negative, 
and the estimates used for ( 3.18| ) are no longer valid. But 



5^i;(|x,-x,|)i^2<5^t;(t,)i^2^ 



(3.43) 



j<i 



LSY 05/08/99 



16 



and because of 2/'(r) + v{r)f{r) > we get 

(3.44) 



So we have the same as in ( p.l8|) . Note that for potentials as in ( |3.42|) / 
satisfies (p.5[), as long as a > 0. 



Remark 3.2 (Homogeneous gas). For the special case of a homogeneous 
Bose gas (i.e. = 0) in a box of volume V, the GP density is simply 

p'^P{x) = N/V = p, (3.45) 

for all X in the box, and the GP energy is given by 

E^^{N,a) =47ia—. (3.46) 

Our method also applies here, if we impose periodic boundary conditions on 
the box. Therefore our upper bound is a generalization of a result by Dyson 
0, who proved an analogous bound for the special case of a homogeneous 
Bose gas with hard-sphere interaction. 

3.2 Upper bound for the chemical potential 

By the same method as in the previous subsection one can derive a bound 
on the increase of the energy when one particle is added to the system. This 
bound will be needed for the derivation of the lower bound to the energy. 

Theorem 3.2 (Upper bound for the chemical potential). Let 

E*{N,a) denote the infimum of the functional 

= y (|V$(x)|2 + \/(x)|<l>(x)|2 + 87ra||$||^|<l>(x)p)rf^a; (3.47) 

with / = A^. Let $* he the positive minimizer of £* (its existence is 
guaranteed by the same arguments as for the GP functional itself), and p* = 
Jq)*^/N. Then 

E^^(N + 1, a) < E'^^iN, a) + E*{1, Na) (l + 0{ap*'^^)) . (3.48) 
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Proof. Let ^I/q^'' be the ground state wave function of H''^\ As test wave 





function for we take 



^(xi, . . . , xn+i) = '^f\xi, Xiv)$*(x7v+i)/(tiv+i), (3.49) 

where / and t^+i are defined as in (|3.30| ) and (|3.4|), i.e., r/(r) is essentially 
the zero energy scattering solution and ^Ar+i is the distance of xn+i from its 
nearest neighbor. We have 

N N 



\i=i i<j 



(3.50) 



For / one uses the estimates 



TV 



/ < 1, /(tjv+i)' > 1 - 5Z (1 - f{\xN+i - x,\f) , (3.51) 



i=l 

and for the derivatives one has 

N 



|VWP = m+i)' = 5^|V./p, (3.52) 



1=1 

and 

AT 



/'(tiv+i)'<5^/'(k^+i-x.|)2. (3.53) 



i=l 
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After division by the norm of \1/ ( p.50| ) becomes 



/|^i^^|2(Ar_iV$*2 J(i_/2)) 

AT 



+ 



/|<)p(iV-iV||$*||^/(l-/2)) 

(3.54) 



^'o^^ does not depend on xat+i. One integrates first over xat+i and tlien over 
the remaining variables. In analogy with the estimates ( |3.33D and ( |3.39D one 
gets 



(2/'(|x|)2 + v{\x\)f{\x\f)d''x < Stto (1 + 0{ap*^/^)) , (3.55) 
and 

WWlo /(I - fi\x\r)d'x < 0{ap*'/'). (3.56) 

This implies 

EQ^(iV + 1, a) < EQ^(iV, a) + (1 + 0{ap*'/'')) . (3.57) 

By scaling, E*(N, a) = NE^^l, Na) and (|3^ follows. □ 



Remark 3.3 (Box with Neumann conditions). Eq. (|3.48| ) also holds 



for a box with Neumann boundary conditions. To see this we note that Neu- 
mann conditions give the lowest energy for the quadratic form (\I'|i^*^^^^'*\E'), 
therefore it is possible to use ( p.49| ) as a test function, even if / does not 



fulfill Neumann conditions. If "^[^^ and $* do, the calculation above is still 
valid, since for 

/ |v(^?/)r = - / fg'^'g + / /I v/r (3.58) 

only boundary conditions for g are needed. 

Note also that in the homogeneous case, i.e. ^ = in the box, E*{N, a) = 
2E^P{N,a). 
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4 Lower Bounds 

4.1 The homogeneous case 

In 1^ the following lower bound was established for the ground state energy, 
^hom^ of a Bose gas in a box of side length L with Neumann boundary 
conditions and v of finite range: 

E^^'^iN, L) > 47ra— (1 - CY^^^^) (4.1) 

with Y = a^N/L^ and C a constant. The estimate holds for all Y small 
enough and L/a ^ y-6/i7 (^^ote that this implies ^ F"^/^''). In the 
thermodynamic limit the constant C can be taken to be C = 8.9, but this 
value is only of academic interest, because the error term —CY^^^'^ is not 
believed to reflect the true state of affairs. Presumably, it does not even have 
the correct sign. 



The restriction of a finite range can be relaxed. In fact, ( [4. 1|) holds (with 



a different constant C) for all positive, spherically symmetric v with 

v{r) < const. r~'''^^5+'^) for r large enough, e > 0. (4-2) 

More generally, if 

v{r) < const, r"'-^'*'^-' for r large enough, e > (4-3) 

then ([4.1|) holds at least with the exponent 1/17 replaced by an exponent 
0(e). We prove these assertions in Appendix B. 

In the next section we shall stick to the estimate ( |4.1| ) for simplicity, but 
in the limit iV — > cxd the explicit form of the error term is not significant so 



a decrease of the potential as in (^.3| ) is sufficient for the limit Theorems |0 
and lOl 



4.2 The lower bound in the inhomogeneous case 

Our generalization of (|4. 1| ) to the inhomogeneous case is as follows: 

Theorem 4.1 (Lower bound for the QM energy). Let v be positive, 
spherically symmetric and decrease at infinity like ( |4.2| ). Its scattering length 
is a = tti/N with ai fixed, as explained in the Introduction. Then as N ^ oo 

E'^^iN, a) > ^^P(iV, a) (l - const. N'^/^^) (4.4) 
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for all R large enough, where Ejf^ is the GP energy in a cube with side length 
2R, center at the origin, and Neumann boundary conditions; the constant in 
( |4.4| ) depends only on ai and R. 

Proof. As in the lower bound will be obtained by dividing space into 
cubic boxes with Neumann conditions at the boundary, which only lowers 
the energy. Moreover, interactions among particles in different boxes are 
dropped. Since v > 0, this, too, lowers the energy. For the lower bound one 
has to estimate the energy for a definite particle number in each box and 
then to optimize over all distributions of the particles among the boxes. 

Step 1 (Finite box): The first step is to show that all the particles can be 
assumed to be in some large, but finite box. Since 



tends monotonically to oo with i?, one knows that the energy of a particle 
outside a cube Kr of side length 2R and center at the origin is at least K{R). 
Hence 



where E^^^^N — n, a) denotes the energy of N — n particles in A/j, with 
Neumann conditions at the boundary. We now apply Theorem p.2| (which 
holds also in a cube with Neumann conditions). Applying the theorem n 
times and noting that E*{1, Na) is monotone in A^ we have 



Ef{N -n,a)> Ef'iN, a) - nE*{l, Na) (l + 0{ap*'/^)) . (4.7) 



Hence there is a constant K' (that depends only on Na), such that K{R) > 
K' implies that the infimum is obtained at n = 0. This is fulfilled for all 
sufficiently large R, independently of A^ if Na is fixed. So we can restrict 
ourselves to estimating the energy in A/j with Neumann boundary conditions. 

Step 2 (Trading V for —p^): We shall now use the GP equation to 
eliminate V from the problem, effectively replacing it by — Svrap^^. We 
write the wave function in A^ as 



K{R) = inf V{x) 



x\>R 



(4.5) 



EQ^(Ar, a) > inf {Ef{N - n, a) + nK{R)}, 



(4.6) 



TV 



...,xn) = f{xi, ...,xn)YI "^fi^i)^ 



(4.8) 



i=l 



LSY 05/08/99 



21 



where denotes the the minimizer of the GP functional in A^j; since it is 
strictly positive, every wave function can be written in this form. Note also 
that and / obey Neumann conditions. We have 

TV 

1=1 k=l i<j 

(4.9) 

where the integrals are over A;^. Using the GP equation ( p.4|) this becomes 

N ^ N / / i-1 



i=l k=l \ \ j=l J 

(4.10) 

Inserting the value ( pl5|) for fiji gives 

- E'i^ = ATrapj^N + Q{f) (4.11) 



Qif) = E 



N / i-1 

J, I U Pfi^k) \VJV + E v{\x, - x,|)|/p - Srrapfixm' 

AN k=l \ j = l 



N 



(4.12) 

Step 3 (Division into boxes): Q{f) is a normalized quadratic form on 
the weighted L^-space L'^{-^R,Y[k=i pW and can be minimized by 
dividing the cube into smaller cubes with side length L, labelled by a, 
distributing the particles among the boxes and optimizing over all distri- 
butions. We therefore have 

infQ(/)> inf VinfQa(/„), (4.13) 

/ {ria} fa 
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where the infimum is taken over all distributions of the particles with 
Y^na = N, and Qa{f) is given by 

a k=l \ j=l 



QM) = Yl 

a k=l 

(4.14) 

where the integrals are over in the box a, k = 1, . . . , 11^. Note that here 
/ = f{xi,..., Xn^), and ( |4.14| ) is the same as ([4.12|) with replaced by Ua 
and Aji with the box a. 

We now want to use ( [4 .11 ) and therefore we must approximate p^^ by 
constants in each box. Let Pa.max and Pa,min, respectively, denote the maximal 
and minimal values of p'^ in box a. With 

...,XnJ = f{x,,...,XnJll ^f{xk), (4.15) 



k=l 



one has 



JTUWi^' (4 16) 



pa, max J |^ | 

This holds for all i, and if we use p^{xi) < Pa,max in ( |4.14[ ), we get 

Qaif) > ^^E'-^^in^, L) - 87rap„,^^n«, (4.17) 



Pa., 



max 



where E^°™ is the energy in a box without an external potential. 

Remark: If we had not taken Step 2 and used instead the division into 
boxes directly on the original Hamiltonian ( |1 . 1| ) we would be considering the 
minimization of 

^E'^°'"(n„,L)+K,^i,n,. (4.18) 
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Such a procedure, however, would not lead to the GP energy. To see this, 
consider the special case of no interaction, i.e., v = and hence also also 
E^"^{na, L) = 0. The minimum of ( [4. 181 ) is then simply N minr^ V {x) , 
whereas the GP energy is in this case times the ground state energy of 
-V^ + V. 

Step 4 (Minimizing in each box): Dropping the subsidiary condition 
'Y^Ua = N can only lower the infimum. Hence it is sufficient to minimize 
each Qa separately. To justify the use of ( ^Af ), we have to ensure that is 
large enough. But if the minimum is taken for some ?7,q,, we have 

(i5;"°'"(n„ + 1, L) - E'-^^in^, L)) > Svrap^.^ax, (4.19) 

Pa, max 

and using Theorem p.2| , which states that 

£;^°"^(n, + 1, L) - L) < 8™-^ (l + 0(n„aVL=')) , (4.20) 



we see that is at least ~ Pa,ma,xL^- We shall later choose L ~ A^~^/^°, so the 
conditions needed for (|4.1| ) are fulfilled for large enough, since Pa,raa.x ~ N 
and hence a ~ N^/^^ and Ya ^ N ^. Thus we have (for large 

enough A^) 

Qaif) > 4vra fP^!^ (1 - CY^^^) - 2n.p^,^.^ . (4.21) 

We now use = a^na/L^ < a^N/Li^ = Y, and drop the requirement that 
Ha has to be an integer. The minimum of (|4.21| ) is obtained for 

This gives for ( [4.11| ) 

AnapRN - A-aa ^ p^^^nin^^ ( ^^TT 

^ \ r a,min \ 



(4.23) 



Now p^^ is differentiable by Lemma |A.6| , and strictly positive. Since all the 
boxes are in the fixed cube there are constants C < oo, C" > 0, such 
that 

Pa,max " Pa.min < NC'L, p,,^i„ > NC". (4.24) 
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Therefore we have, for Y and L small, 



Lnr, < 1 + DY^^'' + D'L (4.25) 
with suitable constants D and D'. Also, 

^naJ^PlMnL' < 4™ /" (pgP)^ < i?^P(iV,a), (4.26) 

and hence 

EQ^(iV, a) > E^^{N, a) (l - DF^/i^ - D'L) . (4.27) 

As last step it remains to optimize the length L. Recall that Y = a^N/L^ 

and Na is fixed. The exponents of in both error terms in (|4.27| ) are the 
same for 



L ~ aAT^/io _ ^-1/10^ ^4 28) 

The final result, therefore, is 

E'^^iN, a) > E^^'iN, a) (l - D"N-^/^°) . (4.29) 

□ 

4.3 The limit theorems 

By Theorems and [4.1| we have (with a = ai/N) 

E^P(Ar, a) (1 + 0(Ar-2/3)) > EQ^(iV, a) > E^^{N, a) (l - 0{N-^/^')) . 

(4.30) 

Dividing by A^ and taking the limit A^ ^ oo we have 

E«P(1, a,) > hm ^^''(^'"i/^) > Efil, a,) (4.31) 

A"— >oo iV 



for all R large enough. Using Lemma 2.2 and taking the limit i? — > oo we 



finally prove Theorem 1.1. 



The convergence of the energy implies also the convergence of the densi- 
ties: We replace V hy V + 5W with W G and denote the corresponding 
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energies by Es{N, a). It is no restriction to assume that l^ + ^W^ > for small 
\S\. Ef^{N, ai/N)/N is concave in S (it is an infimum over linear functions), 
and converges for each 6 to Ef^ll, ai) as — > oo. This implies convergence 
of the derivatives and we have (Feynman-Hellmann principle) 



^ F 



QM 



(AT, a) 



<5=0 



WpTa, 



d_ 

m 



Ef^iN,a] 



5=0 



with given by (|1.3D. In the weak sense we thus have 



(4.32) 



PSa,(a;), 



(4.33) 



which proves Theorem 1.2 



5 Conclusions 

We have proved that the GP energy functional correctly describes the energy 
and particle density of a Bose gas in a trap to leading order in the small 
parameter pa^ (where p is the mean density and a is the scattering length) 
in the limit where the particle number N tends to infinity, but a tends to 
zero with Na fixed. 



Appendix A 



In this appendix we prove Theorem The proof is split into several lem- 
mas. 

Lemma A.l (Strict convexity). For p > 0, ^/p E V, S'^^\^] is strictly 
convex in p. 

Proof. The second term in ( [1.4|) is linear, the third quadratic in p. So it 
suffices to show that the first term is convex. Let pi and p2 be given, with 
$1 = p]^"^ and $2 = P^"^ in T^n- Then also $ = {api + (1 — a)p2)"'^^^ E 
for all < a < 1. We have 

$V$ = a$iV$i + (1 - a)$2V$2 

= (ai/2$i) (a^/^V^i) + ((1 - aY/'<i>2) ((1 - «)'/'V<l>2) 
< + (1 - a)$^)'^' («|V$i|' + (1 - a)|V$2|')'^' 
= $ (a|V$i|^ + (1 - a)|V<l>2n^^^ • 
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Hence 



V$P < a|V<l>iP + (1 - a)|V$2'^ 



□ 

Remark. Because ^ > 0, is also convex in $ G P. But since the 

domain "Dtv is not convex, it is necessary to consider p \—>- S'^^[^/p]. 

Lemma A. 2 (Minimizer). For all N there exists a minimizing $oo ^ T^n, 
with S'^^l^oo] = -E'^^(iV, a). Moreover, |$ooP is unique. 

Proof. Let be a minimizing sequence in T>m, i.e., lim„^oo ^^^^['^'n] = E'^^ . 
It is clear that there exists a constant C, such that ||V$ri||2 < C, H^nlU < 
C and J|$„pl^ < C for all n (recall that V is nonnegative) . Hence the 
sequence belongs to a weakly compact set in L^, as well as in the Sobolev 
space if^ = {$ : ||$||2 + ||V$||2 < oo}, and in the space Ly, defined by the 
norm ||$||y = (/ |$(x) Thus, there exists a $oo e P and a 
weakly convergent subsequence, again denoted by $„, such that 

^ $oo in n n Ll 
V$„ ^ V$oo in L\ 

Because the norm, the Sobolev norm, and the Ly norm are all weakly 
lower semicontinuous, we have 

liminf£GP[$„] >^«P['^'oo], 



|2 



oo 



and it remains only to show that $oo &T^n- Since converges to |$ 
in L\^^ it is clear that H^ooHi ^ Moreover, 

[ |$ooP < ll^ooll^ 
IB J B 

for all bounded regions B. If H^ooHi = N — e with e > 0, then there exists a 
constant Mb for all 5, such that 

for all n > Mb- Since lim|^|^oo V{x) = oo, this would imply J ^ cxo, 

which is impossible because is a minimizing sequence for the functional 
Hence ||$oo||i = A^. 
The uniqueness of |$ooP follows immediately from strict convexity. 
Lemma |A.1[ □ 
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Lemma A. 3 (GP equation). Every minimizing $oo satisfies the Gross- 
Pitaevskii equation ( p.4|) . Conversely, every solution to ( p.4D , with fi given 
by (|2.5[), is a minimizer for . 



Proof. Pick a function / G . The stationarity of at $oo implies 
9 



[£^^[^^ + ef]+^^\\^^ + efr^ 



= 

e=0 



with a Lagrange parameter /i to take account of the subsidiary condition 
II $11 2 = A^. With / real valued one obtains 

-V^Re$oo + VRe^oo + 87ra|$oopR-e $oo = yuRe$oo 

and an analogous equation for Im with / purely imaginary. The value of /j, 
is obtained by multiplying the GP equation with $oo and integrating. By the 
same argument, every solution $ to the GP equation satisfies = E^^ 

and is thus a minimizer. □ 



Lemma A. 4 (Uniqueness). The minimizing $00 is unique up to a con- 
stant phase factor. This factor can be chosen so that $00 is strictly positive. 

Proof. Since £^*^^[|$|] < (by an analogous computation as in the proof 

of Lemma |A.1| ), we know that |$oo| is a minimizer and hence a solution to the 
GP equation. It is thus an eigenstate of the Hamiltonian H = —'V'^ + W with 
W = V^ + 87ra|$ooP (recall that |$ooP is unique), and since it is nonnegative, 
it must be a ground state. Since $00 solves the same equation it is also a 
ground state. Now W G Li^^ and lim|a.|_^oo W{x) = 00, so the ground state 
of H is unique up to a phase and without zeros (see 0|, XIII. 47). □ 

The unique strictly positive minimizer is denoted by 

Lemma A. 5 (Exponential fall-off). For all t > there exists an Mt, 
such that $*^^(x) < Mte"*'''!. In particular, G 

Proof. Put W = V + 87ra($^P)^ and let t > 0. The GP equation implies 
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Using the Yukawa potential Yt{x) = (47r|x|) ^exp(— t|a;|) we can rewrite this 

as 

$GP(a;) = - I Y,{x - y){W{y) - - t')<l>'''' {y)d'y. 

Since > 0, and W{y) — fj, — > for \y\ > R with R large enough, we 
also have 



\y\<R 



Now W^^^ e Lf^c, and hence 



|j/l<^? 

□ 

Lemma A. 6 (Regularity). is once continuously differentiable in 

X G M.^, and V$^^ is Holder continuous of order 1. If V ^ C°° , then 
G C°°. Moreover, E'^^{N,a) is continuously differentiable in a and 



hence in N (by Eq. {^), and dE^^{N,a)/dN satisfi 



es 



Proof. The last lemma and the GP equation imply V^$*^^ G L'^^. Thus 
y<|)GP gxists and is Holder continuous (see 0, 10.2). The property fol- 
lows by a bootstrap argument. The differentiability with respect to the pa- 
rameter a may be shown by a Feynman-Hellmann type argument like analo- 
gous statements (e.g. differentiability w.r.t. nuclear charges) in TF theory 
Eq. (|3D follows immediately from (|r|) and ^^^(1, Na) = S^^[^f^^]. □ 



Lemmas |A.1| - [A.6| complete the proof of Theorem 



Appendix B 

In this appendix we show that (|4.1| ) holds for nonnegative potentials satisfy- 
ing (|4.2|), and that a similar estimate with 1/17 replaced by 0{e) holds under 



the condition (4.3). 
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We cut the potential at a finite radius R which, because of f > 0, can 
only decrease the energy. We thus define 

v{r) =v{r)Q{R-r) (B.l) 

and denote the corresponding scattering length hj a < a. Let u be the zero 
energy scattering solution for the potential v (cf. (|1.2|) ) and put 

h(r)=r-^. (B.2) 
^ ^ u'{r) ^ ' 

The difference a — a can be estimated as follows. Since v{r) and v{r) agree for 
r < R, the same holds for the corresponding scattering solutions. Moreover, 
a = h{R). Hence 



oo 



/ h'f ^^ / u{r)u"{r) 
a — a = j h [rjdr = / -— - — ar 



R R 

oo 



(B.3) 



I XL 2? l / / \27 



u{r 

R R 



where convexity of u has been used to derive the inequality. We remark that 
for i? ^ this simple estimate gives the Spruch- Rosenberg inequality |T^ 



a< \ I v{rydr. (B.4) 

Assuming (|4.2| ) one obtains 







a 



1 



a > a I 1 - const. ( ^ ) | . (B.5) 

.Ry 



Eq. ( [4.1|) holds in any case with a replaced by a, and if we we take 
R oc aY^^^^"^^^' with e' > then the error in ( |B.5|) is of higher order than 
the leading error term in ( [4.11) . We have thus established (|4.1| ) under the 



condition ( |4.2| ). If only the weaker condition ( f4.3| ) holds, then the additional 



error term may be 0{Y^^^-^''^ '^')^) 
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To see the significance of condition ( [4 .31) we also estimate a — a from 



below: 



oo 



oo 



a — a > / u{r)u"{r) > | 




/ 



f (r)(r — a) 



(B.6) 



R 



In order that a is finite the last integral must converge, i.e., a slower decrease 
than is not allowed. 
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